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, Abstract 

<N 

C ■) • We first calculate the binding energy, the pionic and electromagnetic coup- 

ing constants of the lowest lying p-wave heavy baryon doublet A c i, A*^ in 
' the leading order of the heavy quark expansion. Then we calculate the two- 

body decay widths with these couplings and compare our results with other 
approaches. Our results are r(A a -> £ c tt, £ c7 , £* 7 ) = 2.7,0.011,0.001 MeV 
(S| ■ and r(A*j -> S c 7r, £ c7 , £* 7 , A c i 7 ) = 33, 5, 6, 0.014 keV respectively. We find 

■ A c i , A*j — * A c7 is strictly forbidden in the leading order of the heavy quark 

expansion. At the order of 0(l/m c ) their widths are 36,48 keV respectively. 

o 



Q-h! 

Oh! 
CD 



PACS Indices: 14.20.Lq, 13.40.Hq, 13.75.Gx 



I. INTRODUCTION 



Now most of the ground state charm baryons have been found experimentally |lj] . Important progress 
has been made in the search of orbitally excited heavy baryons. The ARGUS E687 || and CLEO 
Q) collaborations have observed a pair of states in the channel A+7r + 7r~, which were interpreted as the 
lowest lying orbitally excited states: A c i(2593) with J p = \~ and A* 1 (2625) with J p = |~. The total 
decay width of the A c i(2593) is 3.6^'g MeV while only an upper limit of < 1.9 MeV has been set for 
A* (2625) up to now [Q. Recently there emerges evidence for the 5*7^ with J p = | , the strange partner 
of the A* 1 (2625). Its width is less than 2.4 MeV. In the near future much more data will be expected. 
We will focus on the strong and electromagnetic decays of the A c i doublet since they are the only well 
established states 

There exist many theoretical discussions on this topic. In || the single pion and two pion strong 
decays and radiative decays of the A c i doublet were discussed within the framework of heavy baryon 
chiral perturbation theory. Due to unknown couplings constants in the chiral Lagrangian, no actual 
decay widths were given. Within the same framework the pionic decay widths were calculated assuming 
the heavy quark effective theory is still valid for the strange quark Q. The coupling constants in the 
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chiral Lagrangian were fixed using the p-wave strange baryon decay widths, which were later used to 
predict the strong decays of the p-wave charm baryons . The two pion width of A c i was estimated 
to be around 2.5 MeV, which was comparable to the total one pion width 3.0 MeV. And the decays of 
was suppressed by more than an order [0]. In j?J the p-wave doublet was treated as the bound state 
of the nucleon and heavy meson. It was found that the decays AdjA^ — > A c 7 were suppressed due 
to the kinematic suppression of the electric dipole moment 0. In Q the constituent quark model was 
employed to study the orbitally excited heavy baryons. Sum rules were derived to constrain the coupling 
constants. The light front quark model, together with underlying SU(2Nf) x 0(3) symmetry for the light 
diquark system, was used to relate and analyse the pionic coupling jj^-jl^]. However, the results have 
strong dependence on the constituent quark mass m q . Varying m q from 220 MeV to 340 MeV, the decay 
widths increase by more than a factor of two |l2| . Within the same framework the electromagnetic decays 
of the p-wave baryons were calculated in [|l3| . In jl4|] both strong and radiative decays were calculated 
using a relativistic three-quark model. After this paper was submitted there appears an interesting paper 
discussing the radiative decays of the ground state heavy baryon multiplets in the framework of heavy 
baryon chiral perturbation theory. In some cases the loop corrections yield sizeable enhancement of the 
deca widths Jl5| . 

It will be helpful to extract these pionic and photonic coupling constants at the quark gluon level 
using QCD Lagrangian. We will treat this problem using QCD sum rules (QSR) |l6), which are successful 
to extract the low-lying hadron masses and couplings. In this approach the nonperturbative effects are 
introduced via various condensates in the vacuum. The light cone QCD sum rule (LCQSR) differs from 
the conventional short-distance QSR in that it is based on the expansion over the twists of the operators. 
The main contribution comes from the lowest twist operators. Matrix elements of nonlocal operators 
sandwiched between a hadronic state and the vacuum define the hadron wave functions. In the present 
case our sum rules involve with the pion and photon wave function. When the LCQSR is used to 
calculate the coupling constant, the double Borel transformation is always invoked so that the excited 
states and the continuum contribution can be subtracted quite cleanly. We have calculated the pionic 
and electromagnetic coupling constants and decay widths of the ground state heavy hadrons jl7]-|l9| and 
possible hybrid heavy mesons J2(| . In this work we extend the same framework to study the strong and 
radiative decays of lowest p-wave heavy baryons, i.e., A c i doublet. 

Our paper is organized as follows: Section || is an introduction. In the next section we derive 
the mass sum rule. The light cone sum rules for the pionic coupling constants are derived in Section 
III. Numerical analysis is presented. In Section [jv| we extend the same framework to analyse the 
electromagnetic processes A c i — * E c 7 etc. In Section |v| we discuss the processes A c i, h* cl — > A c 7 and 
compare our results with other theoretical approaches. The last section is a summary. 
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II. THE MASS SUM RULES FOR THE HEAVY HYBRID MESONS IN HQET 

A. Heavy quark effective theory 

The effective Lagrangian of the HQET, up to order I/tuq, is 

C eS = h v ivDh v + -^—)C + -^—S + 0(l/m 2 Q ), (1) 
IrriQ Zuiq 

where h v (x) is the velocity-dependent field related to the original heavy-quark field Q(x) by 

h v (x) =e im ^ x ^^Q(x) , (2) 

is the heavy hadron velocity. K. is the kinetic operator defined as 

K = K {iD t ) 2 h v , (3) 

where = — (v ■ D)v fJ -, with = 9 M — ig is the gauge-covariant derivative, and S is the 
chromomagnetic operator 

5=| C mag (m Q /fi) h v a^G^hy , (4) 

/ / \ \ 3//3o 

where C mag — I s ^ I , flo = 11 — 2n//3. Note the heavy quark propogator has a simple form in 
coordinate space. 

r°° i + v 

<0\T{h v {x),h v {0)}\0>= dtS(x-vt)—-. (5) 

Jo 2 

B. The interpolating currents 

We introduce the interpolating currents for the relevant heavy baryons: 

VaM = e abc [u aT (x)C l5 d»(x)}h c v (x) , (6) 

»fc+(aO - e abc [u aT (x)C^d b (x)H l5 h c v (x) , (7) 

(x) = e abc [u aT (x)C^u b (x)}Trh c v (x) , (8) 

VA C1 (x) = e abc [u aT (x)C l5 d b (x)H l5 Dlh c v (x) , (9) 

< cl (x) = e abc [u aT (x)C l5 d b (x)]VrDih c v (x) , (10) 



where a, 6, c is the color index, u(x), d(x), h v (x) is the up, down and heavy quark fields, T denotes the 



transpose, C is the charge conjugate matrix, T^ v = —g^ + 9t" = 9^ ~ v ^ v " 's "ft — 7m — vv^, 

and v 11 is the velocity of the heavy hadron. 

The overlap amplititudes of the interpolating currents with the heavy baryons are defined as: 

(0|»7A o |Ac>=/A o ttA o , (11) 
(OltjAjAdH/AdUAci. ( 12 ) 

(01^^) = (13) 
(0|77EjS c )=/ Ea n So , (14) 

(0|^,|S*) = ^<», (15) 

where , are the Rarita-Schwinger spinors in HQET. In the leading order of HQET, /s c = /s* 
and /a c1 = /a*j due to heavy quark symmetry. 



C. The Aqi mass sum rules 

In order to extract the binding energy of the p-wave heavy baryons in the leading order of HQET, 
we consider the correlators 

i J d^e^iOlTirjA^fjA^Om = 1%)^ , (16) 

with u> = k ■ v. 

The dispersion relation for II(w) reads 

n( w ) = / p{s) . ds , (i7) 

J s — u> — te 

where p(s) is the spectral density in the limit ttiq — > oo. 
At the phenomenological side 

P 

n(o;) = — h continuum . (18) 

A c i - uj 

In order to suppress the continuum and higher excited states contribution we make Borel transformation 
with the variable uj to ([IT]) . We have 

n A Acl ["° 

fl cl e-— = / p(s)e--ds , (19) 
Jo 
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a a c1 = rs • ( 2l ) 



where Aa c1 is the A c i binding energy of in the leading order and so is the continuum threshold. Starting 
from so we have modeled the phenomenological spectral density with the parton-like one including both 
the perturbative term and various condensates. 

The spectral density p(s) at the quark level reads, 

where a = -4n 2 (qq) = 0.55GeV 3 , (g 2 s G 2 ) = 0.48GeV 4 , (qg s a ■ Gq) = m 2 (qq), and m 2 = 0.8 GeV 2 . An 
interesting feature of ( pp| ) is that the gluon condensate is of the opposite sign as the leading perturbative 
term, in contrast with the ground state baryon mass sum rules. This may be interpreted as some kind of 
gluon excitation since we are considering p-wave baryons. In the present case the gluon in the covariant 
derivative also contributes to various condensates. 

Two common approaches exist to extract the masses. One is the derivative method. 

L s ° sp(s)e~^ ds 
Io° p(s)e-^ds 

The other is the fitting method, which involves with fitting the left hand side and right hand side of Eq. 
(|l9| ) with the most suitable parameters Aa c1 , f\ c i > s ° m the working region of the Borel parameter. With 
both methods we get consistent results, 

A Acl =(1.1 ±0.2) GeV, 
/ Acl = (0.025 ± 0.005) GeV 4 , 

4 oi = (1.45 ±0.2) GeV (22) 

in the working region 0.5 — 1.3 GeV for the Borel parameter T . For later use we also need the mass and 
overlapping amplitude of the S, A heavy baryon doublet, As c , Aa c , fe e , f\ a in the leading order of a s 

A Ec = (1.0 ±0.1) GeV , 
/ Sc = (0.04 ±0.004) GeV 3 , 

4 C = (1-25 ± 0.15) GeV (23) 

A Ac = (0.8 ±0.1) GeV , 

f Ac = (0.018 ± 0.004) GeV 3 , 

s° A = (1.2 ±0.15) GeV (24) 



III. LCQSR FOR THE PIONIC COUPLINGS 
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A. The correlator for pionic couplings 



We introduce the following amplitudes 

M(A c i -> E c 7r) = g s u^ a u Acl , (25) 

M(A* cl -> £ c tt) = \^ ffd n So7 5^g<. , (26) 
M(A cl -> £*tt) = Viglu^q^qu^ , (27) 

M(A c *j -> £ c *7r) = ^ [<^*„ + 3^(g* ^ - i^JK;, , (28) 

where g = 9^7^, 9^ is the pion momentum. Only the first two decay processes are kinematically allowed. 
Due to heavy quark symmetry, g' s = g 8 , = g" 2 , = in the limit of rriQ — > 00. In other words there are 
two independent coupling constants correpsonding to s-wave and d-wave decays. Note we are unable to 
determine the sign of g s and g d - And we are mainly interested in the decay widths of the p-wave heavy 
baryons. In the following our convention is to let both couplings be positive. 
We consider the following correlators 

d*x e^(0\T( VAcl (x)^M) k(g)) = ±±*G a (u,u') , (29) 



d 4 ' 



x e*-*(Q\T (^(a:)fe e (0)) |tt(q)) = ^g^lT^e^X) , (30) 

where k' = k — q, u> = v ■ k, u>' = v ■ k' and q 2 = m 2 = 0. 

The function G Sj d(u>,u>') has the following pole terms from double dispersion relation. For G s we 

have 

f\ cl fcg s c d 



(A Acl - w')( A s c ~ oj) A Acl -w' A Sc -w 



(31) 



B. Pion light cone wave functions 

To go futher we need the two- and three-particle pion light cone wave functions p3[ : 

< ^( g )|d»7 M75U (0)|0 > = -iUq^ f du e iu " x (<p n (u) + x 2 9l (u) + 0(x 4 )) 

Jo 

+ Ufa - ^) / du e iu «*g 2 (u) , (32) 
1 X Jo 



< n(q)\d(x)ij 5 u(0)\0 > = J7r n / du e mqx <p P {u) , (33) 

m« + m d Jo 



'6(m u + m d ) J 



< 7T(q)\d(x)a^ l5 u(0)\0 > = i( qfl x v - frXn) ^* s I du e luqx Va {u) . (34) 
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< TT(q)\d(x)a a f3j 5 g s G fJ ,„(ux)u(0)\0 >= 

ihAi% ( lc t 9v0 - qvqag^p) - {q^qpgua ~ quqpg^a)] I ^■ i - n {a i )e iqx{ai+va3) , (35) 



< n(q)\d(x)j )ll j 5 g a G al 3(vx)u(0)\0 >= 



+f n -^(q a x f} - / Pa^,, (a^e 1 ^^ 1 ^) (36) 



and 



< 7r(9)|d(a;)7 M 5 s G' ct/ 3(wa;)M(0)|0 >= 

+iU J ^(q*x /3 - qpx a ) f P^ii^je^ ^) . (37) 
q ' % J 

The operator G Q /3 is the dual of G a ^: G Q( a = ^e a( g5 p G 5p ; Po^ is defined as Pa^ = <iaiG?a2<ia3(5(l — 
«i — a2 — as). Due to the choice of the gauge x^A^x) — 0, the path-ordered gauge factor 
Pexp (ig s dux^A^^ux)^ has been omitted. 

The wave function ip^(u) is associated with the leading twist 2 operator, gi{u) and g2 (u) correspond 
to twist 4 operators, and ipp(u) and <Pa(u) to twist 3 ones. The function (p^ is of twist three, while 
all the wave functions appearing in eqs.(p6|), ( |37| ) are of twist four. The wave functions tp(xi,n) (a* 
is the renormalization point) describe the distribution in longitudinal momenta inside the pion, the 
parameters Xi (Y]. ir, = 1) representing the fractions of the longitudinal momentum carried by the quark, 
the antiquark and gluon. 

The wave function normalizations immediately follow from the definitions (|32])-(|3lj): du ^^(u) = 
fydu ip a (u) = 1, f Q du gi(u) = 5 2 /l2, JVaitp±(ai) = $Vanp\\(ai) = 0, fT>ai<px(ai) = 
— J Doii<p\\(ai) = S 2 /3, with the parameter 5 defined by the matrix element: < ir(q)\dg s G aiJj 'y a u\0 >= 

i^Uq^- 

C. The pionic sum rules 

Now the expressions of G s , Gd at the quark level read, 

6/ijr 



G.fauS) = if v J dt J due^-^e^i^eppiu) + (3v>» + [ti^(u)]") 



— (qg s cr ■ Gq))g 2 (u) + ((qq) + j^(q9s(J ■ Gq)) — 2 



4j/a,r / f / du(l-u) / Am'- 1 1 >./•<')- - r) : Vu + «",)],.-.,,(/«;) 



"§/3» / f / ^e^-^^Me^'^^Cg.^WCai) , (38) 



o 
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G d (u, J) = iU I dt [ duue^-^e^i-^tpriu) 



JO 



I,, . i 2 



+ 3«??> + Y6^ 9s(T ' G<z ^ + <2 [ G2 (") + 5i («)])} 

-^/ 3 .y ^ J du(l-u) J Pa,e i - i [ 1 -( ai +" a3 )] e ^'*( Q!1 + ua3 )[l-it(g-t;)(a 1 +na3)]^3 7 r(a i ) 



i 2 Jo 

where ^ = 1.65GeV, / w = 132MeV, F'(u) = and F"(w) = There are two three particle 

terms in the form of (pz-n hi (|3^) > ( |39| ) . The gluon arises from the light quark propagator in the first term 
and from the covariant derivative in the second term. For large euclidean values of u> and uj' this integral 
is dominated by the region of small t, therefore it can be approximated by the first few terms with lowest 
twists. 

After Wick rotations and making double Borel transformation with the variables to and to' the 
single-pole terms in ( |3~i| ) are eliminated. Subtracting the continuum contribution which is modeled by 
the dispersion integral in the region s, s' > sq, we arrive at: 

g s fA cl h c = {$e AA < ASc {6^ P (u )T 5 U(f) + ^(2/ 3 - h ~ h)T 5 U(f) 

+ f (3^(uo) + [mfcr(u)C) T 5 h(f ) + (1 - ^)T 3 / 2 (f ) 

- f (fl2(«o) + |[«Ga(«) + u$i(«)C) (1 ~ i^)T/o(f )} , (40) 

n fc 

where f n (x) = 1 — e _:l: 52 p- is the factor used to subtract the continuum, sq is the continuum threshold. 

k=0 

uq = jr^jj, T = Ti+T 2 ' "^l' ^ 2 are Borel parameters. The functions are defined below. In 
obtaining ( f4(i|) we have used the Borel transformation formula: B^e auJ = S(a — 4) and integration by 
parts to absorb the factors (q ■ v) and l/(q ■ v). In this way we arrive at the simple form after double 
Borel transformation. 
Similarly we have: 

9df Kl f^ = ^e^W^{^ Uo ^( Uo )T 3 /2(f ) - *£{h +h + h)T 3 f 2 (f) 

+ ^w ^(«o)(l ~ i#2 )T/ (f ) - ^u Q [G 2 (u ) + ffi(«o)](l - 1S2)} ■ (41) 

The functions G2(uq), Ii are defined as: 



no 



G 2 (uo) = - / 92(u)du , (42) 





h = / dai / da 2 5 — ^("t) , (43) 

Jo Jo a i 

h = dai I da 2 ^5 -^(a*) , (44) 

Jo Jo a 3 



8 



/ - [ u ° ri d i f^a^ii 1 z ai z a 3' a 3) ii 

-*3 — / dOi\ j [ \\a3—UQ~ai 







da 3 a 3 



"° , (p3n(ati, 1 - u ,u - ai) Z" 1 UQ <^ 37r (u ,a 2 ,l - "o - "2) 

ctai ? -2 +/ da 2 ^ ^ , (45) 

(uo-aii) Jo (1- 1*0-02) 

^4 — / [ j J|ai=li 

Jo a 3 aai 



, <^37r(ai, 1 - Uo,«o - "i) / , ^37r(wo,a2,l-fio-a2) , 

dai ^ / da 2 ^ , (46) 

(uo-air Jo (l-w -a2) 



h = - I u y> 3ir (u , 1 - ito - a 3 ,a 3 ) + / da\ I da 2 5 "^(ai) , (47) 

Jo "3 Jo JO a 3 

j = d[a 1 ip ilT {ai, 1 - ai - a 3 , a 3 )] MQ 

6 dai la3=l-«o 



/•l— «o ^2 a 

/ da 3 - — n[<P37r(ai, 1 - ai - a 3 , a 3 ) — ]|ai=« 
Jo « a i "3 



/ ^"3 / 

Jo Jo 







+ [^37r(ai,l- "1 " «3,a3) 2 JU=u 

a 3 

dai- — [^37r(ai, 1 - ai - a 3 , a 3 ) = — 

aai a 3 

+ [^ 3 7r(ai,l -ai -a3,a3) 3 2 ^ las^T-uo 

a 3 

1 "™ D , d a 3 -ai 

" a 3 j — ¥>37r(ai, 1 — ai - a 3 , a 3 ) 5 — Ui=«o 

da\ ai 



-2f M ml -2 [""das r"*** d r^Kl-Oi-«3,a 3 ) 



-2 jT da 3 ^°"° 3 ^^(ai.l-ai-as.aa) 
+ 2 H** ^ (faa ^(l- a2 -'a3,a 2 ,a3) (48) 

Jo Jo «3 

where a 3 , a.\ are the longitudinal momentum fraction of gluon and down quark inside the pion respectively. 



D. Determination of the parameters for pionic sum rules 

The mass difference between A c i and S c is only about O.lGeV in the leading order of HQET. And 
the values of the Borel parameter T\, T 2 is around 2 GeV in the working region. So we choose to work at 
the symmetric point T\=T 2 = 2T, i.e., uq = which diminishes the uncertainty arising from the pion 
wave functions and enables a rather clean subtraction of the continuum contribution. 

The pion wave functions and their values at the middle point are discussed in ]2^-[25|] . At the scale 
fi = l.OGeV the values of the various functions appearing in (|IJ)-(|d]) at uq — \ are: ^(uq) — (1.5±0.2) 
@, <p P (u ) = 1.142, ^(uq) = 1.463, 9l (u ) = 0.034GeV 2 , G 2 (u ) = 0.02GeV 2 @, ^(«o) = 0, 
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32<uo) = 0, [uipn(u)]'^ Uo = [uLpv{uy^ =Uo = -6, [ugi(u) + uG 2 {u)\l= Uo = -0.29, h = 1.17, h = 1.17, 
h = 31.9, h = -31.9, h = -1-64, h = 247.5, f 3n = 0.0035GeV 2 . We have used the forms in @ 
for Lp^ipti) to calculate integrals I{. The three particle wave functions arc known to next order in the 
conformal spin expansion up to now. The second derivatives need knowledge of the detailed shape of the 
pion wave functions at the middle point. Various sources indicate ^p^iu) is very close to the asymptotic 
form ^{| , which is exactly known. Based on these considerations we have employed the asymptotic forms 
to extract the second derivatives for <p a {u) and ip^iu). 



E. Numerical analysis of pionic sum rules 

Note the spectral density of the sum rule (fi^)-(f4l|) is either proptional to s 2 or s 4 , the continuum 
has to be subtracted carefully. We use sq — (1.3 ± 0.15) GeV, which is the average of the thresholds of 
the A c i and E c mass sum rules. The variation of g s .d with the Borel parameter T and sq is presented 
in Fig. 1 and Fig. 2. The curves correspond to sq = 1.2, 1.3, 1.4GeV from bottom to top respectively. 
Stability develops for these sum rules in the region 0.5 GeV <T<1.5 GeV, we get: 

9sh cl h - (0.5 ± 0.3) x 10- 3 GeV 7 , (49) 
9dhi /e = (2.8 i 0.6) x 10~ 3 GeV 5 , (50) 

where the errors refers to the variations with T and sq in this region. And the central value corresponds 
to T = lGeV and s = 1.3GeV. 
Combining (|§), <M) we get 



g a = (0.5 ± 0.3) , (51) 
g d = (2.8 ± 0.6)GeV~ 2 . (52) 

We collect the values of the pionic couplings from various approaches TABLE I. Note in our notation 3gd 
corresponds to those in (m) . 

We use the following formulas to calculate the pionic decay widths of p-wave heavy baryons. 

r(A cl ^s c7 r) = |i^M , (53) 

27T mA cl 

r(A: i ^E c ^) = ^^|g| 5 , (54) 

where \q\ is the pion decay momentum. We use the values m\ cl = 2.593 GeV, toa^ = 2.625 GeV, 
ms c = 2.452 GeV 0]. In the A c i decays due to isospin symmetry violations of the pion and S c multiplet 
masses, the pion decay momentum is 17, 23, 32 MeV for the final states X+ + 7r~, S°7r + , S+7r° respectively. 
This effect causes significant difference in the decay widths, which are collected in TABLE II. Summing 
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all the three isospin channels we get I\A c i — ► E c 7r) = 2.7 MeV and T(A* 1 — > E c 7r) = 33 keV. The later 
is nearly suppressed by two oders of magnitude due to d-wave decays. 

/From TABLE II we see that our results are numerically close to those from fixing the unknown 
coupling constants from the p-wave strange baryon strong decay widths assuming heavy quark effective 
theory could be extended to the strange quark case |(|. The values of d-wave decay widths from the 
above approach and ours are much smaller than those from the quark models |TJ|,|l(|[l2j . As for the 
s-wave decays various approaches yield consistent results. 



IV. RADIATIVE DECAYS OF P-WAVE HEAVY BARYONS 
A. The correlator 

The light cone photon wave functions have been used to discuss radiative decay processes in 
2|,|^j28H|Jl| in the framework of QCD sum rules. We extend the same formalism to extract the 



electromagnetic coupling consants for the Aqi doublet decays. 

The radiative coupling constants are defined through the following amplitudes: 

M(A cl -> £ c7 ) = eep vp ^e v *u^ [f s cf t a + f d q a vP]^u Acl , (55) 
M (Ad - S c * 7 ) = VSeep^e"*^. [flg pa + /j^hsT^ , (56) 
M(A cl - £ c7 ) = VSee^q V*u Ec [f s g p t a + / d V</] 75 7X cl » (57) 

M(A* cl - £ c7 ) = 3ee^ w ^e^u|. [f s g p t a + / rf V«<K cl , (58) 

where e^(A) and are the photon polarization vector and momentum respectively, e is the charge unit. 
Due to heavy quark symmetry, we have f} = f% = — f s , = f\ = /J — fd- As in the case of pionic 
couplings there are only two independent coupling constants associated with the El and M2 decays. 
We consider the correlator 



/ 



d*x e- tt -*( 7 (g)|T (r? Acl (0)rte,») |0) = e^^^V*^, u,')g pa + F d (u,,u,')q a v p } , (59) 



F Si d{to,uj') has the same pole structures as G Si d{w, w'). 

The light cone two-particle photon wave functions are [ p6| : 

< 7 (g)|g(a:)cr A(t/ g(0)|0 >= ie q e(qq) [ due mqx {(e^ - e v q^[x(/>(u} + x 2 /ii(u)] 

Jo 

+ {(qx){e l _ l x v - e v x^) + (ex)(x^q u - x^) - x 2 (e ll q u - e„^)]ft 2 (u)} , (60) 



< 7(5)I?(z)7m75<?(0)|0 >= f -e q ee lxuprT e v q p x° \ due mqx ^{u) . ((il ) 
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The 4>(u), ^{u) is associated with the leading twist two photon wave function, while gi{u) and g 2 {u) are 
twist-4 PWFs. All these PWFs are normalized to unity, J Q du f(u) = 1. 

We want to emphasize that the photon light cone wave functions include the complete perturbative 
and non-perturbative electromagnetic interactions for the light quarks in principle. Yet the interaction 
of the photon with the heavy quark is not parametrized and constrained by the photon light cone wave 
functions. It seems possible that the photon couples directly to the heavy quark line. This is different 
from the QCD sum rules for the pionic couplings since pions can not couple directly to the heavy quark. 
However the real photon coupling to heavy quark involves a spin-flip transition, which is suppressed by 
a factor of 1/itiq 0. So it vanishs in the leading order of I/thq expansion. Since we are interested in 
the leading order couplings f s<d , it's enough to keep the photon light cone wave functions for the light 
quarks only. 

Expressing ( |59| ) with the photon wave functions, we arrive at: 

e d )(qq) J dt J due^^e^ii-x^u) 

+ - h 2 (u))] + ^/V(u)t} + • • • . (62) 



F s (ui,u>') = \ {e u 



/■oo />1 -i 

F d {u>,u>') = \{e u - e d )(qq) / dt / duS^^e^'u^-x^u) 
7T ,/n .In t 



+~(hi(u) - h 2 {u))\ + ^Mu)t} + ■■■ . (63) 



The final sum rules are: 



fsfA cl fc c = -£j( e « - e d )e- k ^ Jk {x0(«o)r 5 / 4 (f ) 

-Muo) - M^o)]T 3 / 2 (f ) + ^Mua^fo^)} , (64) 



Ob A -y ' S c "^0 

/d/Aci/Ec = -j-f (e« - e d )e a ^ c Mo { x ^( Mo )T 4 / 3 (— ) 

-[hi(u ) ft 2 (u )]T 2 /i(f ) + ^M«o)} • (65) 



B. Numerical analysis of the photonic sum rules 

The leading photon wave functions receive only small corrections from the higher conformal spins 
so they do not deviate much from the asymptotic form. We shall use f27|| 

<f>(u) = 6u(l - u) , (66) 
= 1 , (67) 
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fti(«) = ~(l-«)(3-u), (68) 

h 2 {u) = -\(\-uf . (69) 

with / = 0.028GcV 2 and X = -4.4GeV 2 ||-fi at the scale p = lGeV. 

The variation of f s ,d with the Borel parameter T and so is presented in FIG. 3 and FIG. 4. Stability 
develops for the sum rules (|64|), ( |65| ) in the region 0.5 GeV <T<1.5 GeV, we get: 

/ s /a c1 /s = (2.0 ± 0.8) x 10- 4 GeV 6 , (70) 

fdfA cl h = (4.8 ± 1.2) x 10- 4 GeV 5 , (71) 

where the errors refers to the variations with T and sq in this region. And the central value corresponds 
to T = l.OGeV and s = 1.3GeV. Our final result is 

f s = (0.20 ± 0.08)GeV _1 , (72) 

f d = (0.48 ± 0.12)GeV~ 2 . (73) 

The decay width formulas in the leading order of HQET are 



r(A cl - 


-s c7 ) 


= i6a|ai» mSc [ft - 












r(A cl - 


- Kj) 












r(Ki - 


-s c7 ) 


- 4a|«f ^ [/ 2 + 




cl 




nKi - 




= 20a|#^l[/ s 2 




TO A » 

cl 





(74) 

where = 134, 72, 164, 103 MeV is the photon decay momentum for the above four processes. The d- 
wave decay is negligible. The decay width values are collected in TABLE III. The uncertainty is typically 
about 50%. 

The decays Aqi — > Sq7 do not occur in the leading order in the bound state picture 0. Due to the 
unknown coupling constant crs in the chiral lagrangian for the heavy quark electromagnetic interactions, 
no numerical values are available || . However the decay width ratios of the four final states are exactly 
the same as ours if we ignore the isospin violations of the heavy multiplet masses in the heavy quark 
limit. Our results are much smaller than those from various versions of quark models [^3| , ^4[ , which may 
indicate that the l/m c correction is important. 

V. THE PROCESS A C i -> A C 7 ETC 

As can be seen later the radiative decay processes of p-wave A cl doublet to A c is quite different from 
those in the previous section. We present more details here. The possible El decay amplitudes are 
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M(A c i -> A c 7) = e/i p e*MA c [fff u • g - u'Vfr/iTsUAci , (75) 

M(A* cl -> A c7 ) = v^e/ 1 ; e ;u Ac ' 9 - w'VK^ . (76) 

Due to heavy quark symmetry h p — h' p . 
We consider the correlator 

H = *y" d*a:e <fc -*(7(g)|T (^(1)^(0)) |0) = i±H e J[^t, • g - w'Vh^/y^u/) . (77) 

We first calculate the part solely involved with the light quark, which can be expressed with the 
photon wave functions. We get 

n = 2i J dtj d 4 x e lk x D t 5{x - vt) l5 ^-{Tr[ l5 CiSl{x)C l5 < j(q)\d(x)d(Q)\0 >] + (u «-» d)} , (78) 

where summation over color has been performed. There are two types of terms with even 7 matrices in 
the trace. The first one is connected with ip(u) and the trace looks like Tr [75 Cx T (^757^ 75]. The second 
is involved with (f)(u), hi(u), h,2(u) and the trace looks like Trl-f^ClC^a^]. In both Aqi, Aq states the 
up and down quarks are in the + state, which leads to the presence of 75C and C75 in both traces. 
Clearly both traces vanish. This property results from the underlying flavor and spin structure of the 
light quark sector. In other words the light quark contribution is zero to all orders of the heavy quark 
expansion in the framework of LCQSR with the commonly used interpolating currents (^) and (^|) for Aq 
and Aqi respectively. The decays Aqi — > Aq7 and Aq 1 — > Aq7 happens only when the photon couples 
directly to the heavy quark line. 

Now let's move to the part involved with the heavy quark. At first sight there are two types of terms 
in the leading order of heavy quark expansion. The first one comes from the insertion of the operator 
i J \h v (y)iv ■ Dh v (y)]d A y in (fT7|), which contributes a factor v ■ e*(A) to the decay amplitude. For the real 
photon v ■ e*(A) = so it drops out. The other possible term arises from the covariant derivative in ri\ cl , 
which leads to a nonzero correlator. For the tensor structure ie*^^r- we have 



n(w,t/) = 4 / dte luJt 



. 6 < g 2 s G 2 > a 2 , 

where the photon field has contributed a factor e~ lq ' x . It's important to note only the variable u>' 
appears in (|79|). It's a single pole term which must vanish after we make double Borel transformation to 
the variables to, to' simultaneously. We have shown there is no leading order El transition in ([75]) arising 
from the photon couplings to the heavy quark line in the leading order of heavy quark expansion. Based 
on the same spin and flavor consideration we know that radiative decay processes like Sq! — > Aqx'J, 
Aqi — * Aqi7, Eqi — > Sqi7 are also forbidden in the leading order of 1/toq expansion, where we have 
used notations in ||. 

We may rewrite the decay amplitudes as 
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M(A cl -> A c7 ) = ef p F^u Ac a^ l5 u Aal , (80) 

M(A* cl -» A c7 ) = 2s/3efiF^u Ac u Ki , (81) 

M (A* x -> A cl7 ) = 2%/3e/ p 2 F Ml/ u Acl 7f 75 u A , i . (82) 
Due to heavy quark symmetry we have 

U = fp = fp ■ (83) 

Note f p =\hp. 

In these decays we know the light quarks do not contribute. However, the J p of the light diquark 
changes from l~ to + which ensures the decay A cl — > A c7 is an El transition. The angular momentum 
and parity J p = i + of the heavy quark does not change so the coupling constant f p is the same as that 
for the heavy quark Ml transition, which is induced by the magnetic moment operator 

< 84 » 

Another approach is to consider the three point correlation function for the tensor structure et 7 5 

i J d*xd i ze ik *- ik '*(0\T{T, Acl (x), m ^ (0) ,fy Acl (z)}|0) = U 3 (lj,l;')±±1 , (85) 
with w = k ■ v, u>' = k' • v. 

n 3(^') - £4 f M^+^t^ + , (86) 







After the double Borel transformation and continuum subtraction we get the sum rule for h p 

h p (A Acl - A A Jf Ac J Ac e-^^ = -l^ { 36T 8 /7 (f ) + ^f^T^fA} . (87) 
Dividing @ by © we get 



hp = — ( 



3(A Acl - AaJ/a T 8/ 7 (f ) - <^>TV 3 (f ) + ' 
Numerically we have /i p w which is consistent with (|84|). The decay widths formulas are 



r(A cl ^A cl ) = e 2 c a\Q\ 3 ™ A 



m Aal m 2 c 



T(A* cl - A c7 ) = e^|?f 



2 



m A * m c 

cl <- 

r(A^ - A cl7 ) = e^gf . (89) 

mA* ti^ 

cl 

The decay momentum is 290, 320, 32 MeV respectively. We take m c = 1.4 GeV. The numerical values are 
collected in TABLE III. These widths comes solely from the C(1/toq) correction. But their numerical 



15 



values are greater than those leading order widths for the channels £ c 7, £*7- The reason is purely 
kincmatical. The decay momentum for the final state A c 7 is three times larger. For the p-wave decay 
there appears an enhancement factor of 27. 

These widths in ( |S9| ) are propotional to Therefore the corresponding radiative decays A^i — > A&7, 

— > Afc7, Aj^ — > Ahi7 are further suppressed by a factor (f 4 -^-) 2 ~ 40. The widths of the first two 
decays are around 1 keV. 

If we use naive dimensional analysis to let crt @ hi TABLE III be of the order of unity or simply 
assume that the El transition coupling constant h p in ( |75| ) is of the same order of Ml transition one ||, 
we would get a width 0(100) keV. Our result is in strong contrast with those from the bound state picture 
[§, where r(A c i,A* 1 -> A c7 ) = 16,21 keV and r(A b i,A* 1 -> A h 7) = 90,119 keV. Future experiments 
should be able to judge which mechanism is correct. 

It was noted that the radiative decays Aqi — * Aq was forbidden in the leading order of heavy quark 
symmetry assuming one-body transition operators, which arises from a complete cancellation due to the 
specific spins of light constituent quarks in the antisymmetric initial and final state [ fl3|] . The point is 
consistent with our observation of the vanishing contribution of the light quark sector to this radiative 
process. 

/From our calculation we know the d-wave single pion width of A*j is 33 keV and the estimate in || 
yielded 35 keV for the two pion decay width. It's interesting to notice that the radiative decay widths are 
48,5,6,0.014 keV for the final states A c 7, £ c 7, £*7, A c i7 respectively. The width of the decay channel 
A c 7 is bigger than either of that of the strong decay modes. The A* x should be a narrow state with a 
total width about 130 keV. 

The two pion width of A c i is about 2.5 MeV ||. /From TABLE II and III the one pion and 
electromagnetic widths are L(E c tt, A c 7, £ c 7, £*7) = 2.7, 0.048, 0.011, 0.001 MeV. Its total width is about 
5.4 MeV. 

It's believed that A^i lies below £{,7r, A^mr threshold. If so its dominant decays are electromagnetic. 
From our calculation we see r(Ahi — > A&7, £;,7, ££7) = 1,11,1 keV if we assume the same decay momen- 
tum as in the A c i decays. Its total width is about 13 keV. It will be a very narrow state. Clearly the 
radiative channels £^7 will be very useful to find them experimentally. 

The major decay modes of A^ might be d-wave one pion decay and electromagnetic decays to £;, 
doublet if the two pion mode is not allowed. Their widths are r(A^ — > £;,7r, A&7, £(,7, £j*7) = 33, 1, 5, 6 
keV if we assume the same decay momentum as in the A*j case. It's also a very narrow state with a 
width of 45 keV. 

Before ending this section we want to improve our previous calculation of radiative decays of excited 
heavy mesons jL9]|. (1) First the s-wave terms involved with g s should not appear in (l + ,2 + ) — > (Q~, 1 — )7 
processes. All decays are M2 transitions. The gl in the decay width formulas should be replaced by 
i^|g| 4 . The last eight widths in Eq. (94) should read 2,8,3, 11,6,23,7,27 keV respectively, which is 
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much smaller than original wrong ones. (2) The El transition (0 + , 1 + ) — > (0 _ , 1~)7 decays was identified 
as s-wave decays. This was misleading. The factor (q ■ v) should be in the tensor structure to ensure the 
El transition structure in Eq. (47) in p9[ . We present the correct sum rules for g\ below. 

a A -,i/2 +A +.i/2 s n x 

ffi/-,i/a/+,i/a = - { X 4>(u )T f (^) - 9l (u )-} , (90) 

where s = cj c /2 = (1.5 ± 0.2) GeV. Numerically we have g x = (1.6 ± 0.2) GeV" 1 . 



VI. DISCUSSIONS 

In our calculation only the errors due to the variations of T and sq are included in the final results for 
5s. di fs,d- The various input parameters like quark condensate, gluon condensate, x, / etc a lso have some 
uncertainty. Among these the values of the pion and photon wave functions introduce largest uncertainty. 
Although their values are constrained by either experimental data or other QCD sum rule analysis, they 
may still lead to ~ 25% unceritainty. Keeping the light cone wave functions up to twist four also leads to 
some errors. However the light cone sum rules are dominated by the lowest twist wave functions. Take 
the sum rule ( |65| ) for /<j for an example. At T = 1 GeV, the twist-four term involved with hi, I12 is only 
9% of the leading twist term after the continuum subtraction. In other words the light cone expansion 
converges quickly. So we expect the contribution of higher twist terms to be small. There arc other 
sources of uncertainty which is difficult to estimate. One is the QCD radiative correction, which is not 
small in both the mass sum rule and LCQSRs for the pionic coupling constants of the ground state heavy 
hadrons in HQET. But their ratio depends only weakly on these corrections because of large cancellation 
p5| . Numerically the radiative corrections are around 10% of the tree level result. 

Another possible source is the I/tuq correction for the charmed p-wave baryons. The leading order 
coupling constants g s ^d etc will be corrected by terms like g' s d/ m Q> which will affect decay widths. For the 
charmed hadrons I/tuq corrections are sizable and may reach 30% while such corrections are generally 



less than 10% of the leading order term for the bottom system |17|. Especially for the El transition 
coupling constant f Sl the correction is of the order -rS--, which may be comparable with the leading order 
one for the charm system. One is justified to use these coupling constants to calculate the decay widths 
of the p-wave bottom baryons. Unfortunately data is still not available for the p-wave bottom baryons. 
So we have calculated the p-wave A c i doublet decay widths with some reservation. 

In summary we have calculated the pionic and electromagnetic coupling constants and decay widths 
of the lowest p-wave heavy baryon doublet. We compare our calculation with different approaches in 
literature. We hope these results will be useful in the future experimental search of A51, A^ baryons. 
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TABLE I. Pionic coupling constants 



Coupling 


Our 


Ref. [ 


U 


1 


Ref. Jo 




9s 
3#d 


0.5 ±0.3 
(8.4 ± 1.8)GeV" 2 


0.52 
21.5 GeV" 2 


0.665±0.135 
50.85±14.25 GeV" 2 



TABLE II. Single pion decay widths 





Our 


Ref. § 


Ref. @ 


Ref. @ 


Experiment | 


S-wave transitions 




Aci ; s(2593) -► E°7r+ 
Aci ; s(2593) -> E+tt 
A c i ; s(2593) -> E++7T- 


0.86 MeV 
1.2 MeV 
0.64 MeV 


0.89 ± 0.86 MeV 
1.7 ±0.49 MeV 
0.55±J;! 5 MeV 


0.83 ± 0.09 MeV 
0.98 ±0.12 MeV 
0.79 ± 0.09 MeV 


1.775 ± 0.695 MeV 
1.18 ±0.46 MeV 
1.47 ± 0.57 MeV 


0.86±° 7 5 3 6 MeV 
r(A c i ;S ) = MeV 
0.86l°;^ MeV 


D-wave transitions 




AJ 1;S (2625) -» E°7r+ 
A* 1;S (2625) - E+tt 
A: 1;S (2625) -» E++7T- 


0.011 MeV 
0.011 MeV 
O.OllMeV 


0.013 MeV 
0.013 MeV 
0.013 MeV 


0.080 ± 0.009 MeV 
0.095 ± 0.012 MeV 
0.076 ± 0.009 MeV 


0.465 ± 0.245 MeV 
0.42 ± 0.22 MeV 
0.44 ± 0.23 MeV 


< 0.13 MeV 
r(A^) < 1.9MeV 

< 0.15 MeV 



TABLE III. Radiative decay widths 





Our 


Ref. @ 


Ref. @ 


Others 


Experiment 1 


A cl;S (2593) -> A+7 


0.036 MeV 





0.115 ±0.001 MeV 
0.016 MeV @ 


0.1914t MeV § 


< 2.mto ll MeV 


A cl;S (2593) -> E+7 


0.011 MeV 


0.071 MeV 


0.077 ± 0.001 MeV 


0T274 S [§ 




A cl;S (2593) -» E c *+ 7 


0.001 MeV 


0.011 MeV 


0.006 ± 0.0001 MeV 


0.0064s | 




A; 1 . S (2625) -» A+7 


0.048 MeV 





0.151 ±0.002 MeV 
0.021 MeV @ 


0.2534t MeV §) 


< 1 MeV 


A: 1;S (2625) -» E+7 


0.005 MeV 


0.13 MeV 


0.035 ± 0.0005 MeV 


0.0584s | 




A: 1;S (2625) - E c *+ 7 


0.006 MeV 


0.032 MeV 


0.046 ± 0.0006 MeV 


0.0544s [| 
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Figure Captions 



Fig. 1. Dependence of (/s/a c1 /e c on the Borel parameter T for different values of 
the continuum threshold so. /,From top to bottom the curves correspond to so = 
1.4,1.3,1.2 GeV. 

Fig. 2. Dependence of gdfA'fs c on T, so. 
Fig. 3. Dependence of f a f\ cl fs c on T, s . 
Fig. 4. Dependence of /<j/a c1 /e c on T, s . 
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